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ABSTRACT 

We propose that detailed data on the star formation history of a dwarf spheroidal 
galaxy (dSph) may be used to estimate the evolution of the total mass Mg [t) for cold 
gas in its star-forming disk. Using Fornax as an illustrative example, we estimate its 
Mg{t) and the corresponding net gas flow rate ^F{t) assuming a global star formation 
rate = A*(t)[Mg(t)/MQ]“ consistent with observations of nearby galaxies. We 
identify the onset of the transition in from a net inflow to a net outflow as the 

time fsat at which the Fornax halo became a Milky Way satellite and estimate the 
evolution of its total mass Mh{t) at t < tsat using the median halo growth history 
in the ACDM cosmology and its present mass within the half-light radius derived 
from observations. We examine three different cases of a = 1, 1.5, and 2, and justify 
the corresponding A*(t) by comparing the gas mass fraction fg{t) = Mg{t)/Mh{t) at 
t < ^sat with results from simulations of gas accretion by halos in a reionized universe. 
We find that the Fornax halo grew to Mh(tsa.t) 2 x 10® Mq at fgat ~ 5 or 8 Gyr, 
in broad agreement with previous studies using data on its stellar kinematics and its 
orbital motion. We describe qualitatively the evolution of Fornax as a satellite and 
discuss potential extension of our approach to other dSphs. 

Key words: galaxies: dwarf — galaxies: evolution — galaxies: formation — galaxies: 
individual (Fornax dwarf galaxy) — Local Group. 


1 INTRODUCTION 


As a result of hierarchical structure formation, the cold 
dark matter (CDM) halo associated with a main galaxy 
has many subhalos. In principle, this can account for the 
dwarf spheroidal galaxies (dSphs) orbiting the Milky Way 
(MW). However, due to the complicated gas dynamics in¬ 
volved in star formation and the resulting feedback, it is 
far from straightforward to associate the subhalos found in 
CDM simulations with the actual dSphs possessing stellar 
populations. A prominent example is the so-called too-big- 
to-fail (TBTF) problem: the largest subhalos of an MW-like 
halo appear unsuitable for hosting the brightest MW dSphs 
because of incompatible mass distributions (e.g., Boylan- 


[Kolchin, Bullock fc Kaplinghat|2011[). Proposed solutions to 
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this problem mostly call for either modifying the nature of 
DM (e.g., |Herpich et al.|2014|) or including baryons and the 


associated gas dynamics in CDM simulations (see [Pontzen 
Ifc Governatoll2014l for a review and ISawala et al. 1120141 for 
recent simulation results on the Local Group galaxies). It 
was also shown that there might not be a TBTF problem 
if the CDM halo of the MW has a relatively low mass of 
« 10^^ Mq (e.g., Wang et al. 20121. In any case, gas dy¬ 


namics is the crucial bridge linking DM halos with observed 
galaxies. 

In this paper we propose an empirical approach to study 
halo evolution and global gas dynamics of dSphs. We use 
Fornax, the brightest MW dSph, as an illustrative example. 
Our approach is in the same spirit as that of|Qian fc Wasser^ 


burg (20121, who used the data on metallicity distributions 


of MW dSphs to infer the generic evolution of their gas. We 
show that the global star formation rate (SFR) tp{t) of For¬ 
nax derived from observations by |de Boer et al. ( 2012b| ) may 
be used to estimate its total gas mass Mg(t) and the corre¬ 
sponding net gas flow rate AF{t) as functions of time t. We 
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consider two regimes of gas evolution separated by the time 
tsat, at which Fornax ceased evolving independently and be¬ 
came an MW satellite. We assume = \t{t)[Mg{t)/MQ]°‘ 
with A*(t) = A*(tsat) for t > tsat- We consider that the ac¬ 
cretion of Fornax by the MW caused the former’s global 
gas flow to change rapidly from a net inflow to a net out¬ 
flow and identify the onset of this transition in AF{t) as 
tsat- We assume that the total mass Mh{t) of the Fornax 
halo at f < fsat follows the median halo growth history in 
the ACDM cosmology and determine Mh{tsat) by requiring 
that the corresponding density profile give the mass enclosed 
within the half-light radius as derived from observations. 

Considering the central role of the assumed global star 
formation law (SFL) Tp{t) = A*(t)[Mg(t)/M 0 ]“ in our ap¬ 
proach, we examine three different cases of a = 1, 1.5, and 
2, and justify the corresponding At(t) by comparing the gas 
mass fraction fg{t) = Mg{t)/Mh{t) at f < tsat with the re¬ 
sults from cosmological simulations of gas accretion by halos 
in a reionized universe. For each case we also perform statis¬ 
tical realizations of the data on to gauge how uncertain¬ 

ties in ip{t) affect our results, especially AF{t). Using repre¬ 
sentative results for the baseline case of a = 1.5, we give a 
qualitative description of Fornax’s overall gas evolution, es¬ 
pecially its gas loss through a combination of ram-pressure 
stripping and tidal interaction with the MW at t > taat- 

This paper is organized as follows. In ^we discuss the¬ 
oretical and empirical SFLs and motivate our assumed SFL 
for dSphs. In we use Fornax as an illustrative example 
to estimate Mg{t), AF{t), Fat, and Mh(tsat) from its star 
formation history (SFH), approximately taking into account 
uncertainties in the SFH. We also compare our estimates of 
taat and Mh(taat) with relevant results from previous stud¬ 
ies, justify choices of A*(t) by considering gas evolution at 
t < taat, and examine how our results are affected by our as¬ 
sumed SFLs. In ij^we qualitatively relate AF{t) at t > taat 
to the orbital motion of Fornax as an MW satellite, and 
discuss how it lost all of its gas through a combination of 
ram-pressure stripping and tidal interaction with the MW. 
We summarize our results and discuss potential extension of 
our approach to other dSphs in © 


2 SFLS AND STAR FORMATION IN A DSPH 


An in-depth study of the halo evolution and the associated 
gas dynamics of a dSph would require simulations that incor¬ 
porate hierarchical structure formation in the ACDM cos¬ 
mology and treat a wide range of physical processes, such 
as gas accretion by the hosting halo before and after reion¬ 
ization of the universe, cooling and condensation of gas to 
form stars, conversion of cold into hot gas by radiation and 
supernova (SN) explosion, and gas expulsion from the halo 
(see e.g., Somerville & Dave|2015 for a review). A represen¬ 
tative example of this approach is the recent work of |Shen| 
et al. (20141, who simulated the formation and evolution of 


seven field dwarf galaxies with present-day virial masses of 
Mvir = 4.4 X 10®-3.6 X 10^° Mq. Our approach here is empir¬ 
ical and phenomenological, especially regarding star forma¬ 
tion and gas inflows and outflows. Central to our approach 
is the SFL relating the SFR to gas properties. Therefore, 
we start with a brief discussion of theoretical and empirical 


SFLs (see e.g., Kennicutt fc Evans|2012 Krumholz|2014 for 
comprehensive reviews). 


2.1 Theoretical and empirical SFLs 


Star formation is most directly associated with molecular 


gas (see e.g., Gao & Solomon 2004 and Krumholz & Thomp- 
|son|2007| for discussion of observations using different molec¬ 
ular tracers). For a cloud with a total mass Mmoi of molec¬ 
ular gas (predominantly consisting of H 2 molecules and the 
associated He atoms), a theoretical estimate of the SFR 
(e.g., Krumholz fc Tan|2007 1 is 


dM* 


eff- 


Mniol 


dt ts 

where M* is the mass of stars formed. 


tg = 


Stt 


32G'Pmol 


( 1 ) 


( 2 ) 


is the free-fall time at the density pmoi of the molecular 
gas, G is the gravitational constant, and eg is the fraction 
of the molecular gas turned into stars during one free-fall 


time. 


Krumholz & Tan 


(20071 obtained eg ~ 10 


from 
- 10 - 


observations covering a wide range of rimoi = Pmoi /rrip 
10“^ cm“®, where rup is the proton mass. An alternative form 
of Eq. (Q is 

/H2S9 


Esfr = Cff■ 


tg 


( 3 ) 


where Esfr is the SFR per unit area. Eg is the surface 
mass density of gas in all forms, and /h 2 is the fraction of 
H mass in H 2 molecules. As star formation typically occurs 
in gaseous disks, Eq. © is more convenient to use, espe¬ 
cially for comparison with observations. Both Eqs. 0 and 
^ can be applied to individual star-forming regions and in 
the disk-averaged form, to individual galaxies (see review 
by |Krumholz|20l^ . For illustration, we compare them with 
two observational results below. 

From observations of nearby galaxies on sub-kpc scales, 
|Bigiel et al.| ( |2008| obtained an SFL that applies to local 
regions within a galaxy: 


Mq yr“i kpc 


= 10 ” 


-'H 2 


10 M© pc" 


1.0±0.2 


( 4 ) 


where Eh 2 is the surface mass density of H 2 molecules in 
a star-forming region. The above SFL covers Eh 2 ~ 3- 
50 Mq pc“^. Taking the central values of the fitting pa¬ 
rameters, we can rewrite it as dM^/dt ~ Mmoi/foepi ~ 
1 . 4 MH 2 /fdepi, where Mua is the total mass of H 2 molecules 
and tdepi ~ 2 Gyr is the time needed for star formation 
to deplete all the molecular gas. This depletion time can 
be understood from tdepi ~ tg/eg [see Eq. ([^]. The star¬ 
forming regions studied by |Bigiel et al. ( 2008| were resolved 
to 750 pc. Assuming a reasonable scale-height H ~ 100 pc 
for such regions with a typical Eh 2 ~ 10 Mq pc“^, we ob¬ 
tain Pmoi ~ I. 4 EH 2 /H ~ 0.14Mq pc“® (umoi ~ 5.7 cm“®), 
which corresponds to tg ~ 22 Myr [see Eq. ([^] or tdepi ~ 
2.2 Gyr for eg ~ 10"^ 

Based on observations of 61 normal spiral galaxies and 


36 starburst galaxies, Kennicutt (1998 1 found a tight cor¬ 


relation between Esfr and Eg, both of which are averaged 
over the disk of an individual galaxy. This global SEL is 


© 2015 RAS, MNRAS 000,[THT3l 








































Evolution of gas in Fornax 3 


jSFR 


Mq yr“i kpc 


32 = (2.5±0.7)xl0' 


Mq pC" 


.(5) 


where Eg ranges from ~ 10 to 10"* Mq pc“^ for the galaxies 
used. The power-law form Esfr oc Eg was first suggested 
by [Schmidt] ( |1958| ) , and Eg. (|^ is com monly referred to as 
the Kennicutt-Schmidt law. Kennicutt (19981 also found an 
alternative form of this SFL: 


Esfr 

Mq yr“i kpc“ 


= 1.7 X 10' 


^Mq pC-2^ 


' yr 

V ''"dyn , 


( 6 ) 


where Tdyn is the dynamic timescale taken to be the orbit 
time at half of the outer radius of the star-forming disk and 
ranges from ~ 10^ to ~ 10® yr. Detailed models to explain 
Eqs. it and (j^ were given by jKrumholz fc McKee] ( |2005[ ) . 
Here we consider the limit of Eg > 10^ Mq pc“^, for which 


gas predominantly consists of H 2 molecules (Krumholz & 
|McKee|[2005P and Eg. (]^ can be applied to give Esfr ~ 
efiEg/tff. Evaluating ts at the gas density pg ~ TiglT-L, we 
obtain 

Esfr 


Mq yr“i kpc 


1.2 X 10' 


( 

/100 pc\ 

VlO-V 

\ H ) 


Mq pc- 


(7) 


which is in good agreement with Eg. (J^ for erf ~ 10 ^ and 
disk scale-heights of "H ~ 100 pc. A similar argument was 


used by Kennicutt] ( 19981 to explain qualitatively the power- 


utt|(; 


law form of Eg, 


2.2 Examples of treating star formation in galaxy 
simulations 

One way to treat star formation in modeling galaxies is to 
use Eg. Q or its equivalent, Eg. §. In either case, es¬ 
timates of /h 2 are required. For systems with sufficiently 
high SFRs, /h 2 depends on Eg, the metallicity, and the 


dumpiness of gas (Krumholz, McKee & Tumlinson 20091. 


For systems with very low SFRs, /h 2 also depends on Esfr 
||KrumhoIz]|2013[). Incorporating all the above dependences 


of /h 2 , Krumholz (20131 showed that Eg. (J^ provides a 


good description of the data on star formation for condi¬ 
tions ranging from poor to rich in molecular gas. 

However, many simulations of galaxy formation do not 
evaluate /h 2 in implementin g Eg. ([l| or ([3| for s tar forma¬ 
tion (see e.g., the review by [Somerville fc Davep015| . In¬ 
stead, stars are assumed to form when gas has cooled below 
some temperature Tmax and condensed to some threshold 
numbe r density of atoms rimin. For example, [Stinson et al.[ 
(20071 chose Tmax = 1.5 x 10"^ K and rimm = 0.1 cm“®. For 
a region with a total mass Mg of the gas that satisfies the 
conditions and has a dynamic timescale tdyn oc (Gpg)“^^^, 
stars were assumed to form stochastically at a rate 


^=0.05^ 

(it ^dyn 


( 8 ) 


In comparing the above equation with Eg. 0 , Mg is a proxy 
for Mmol while the numerical coefficient and tdyn play the 
roles of eg and tg, respectively. More recent simulations of 


Shen et al. (20141 adopted Tmax = 10 K, rimin = 10 


and an SFR per unit volume 


dp* 

dt 


= 0.1 


Pg 

tdyn 


(9) 


where p* is the density of stars formed. The above SFL is 
qualitatively the same as Eg. (Jm. 


2.3 Assumed dSph SFL based on gas mass 

As our goal is to gain some understanding of the overall pic¬ 
ture for halo evolution and gas dynamics of dSphs in general 
and Fornax in particular, we are interested in the relation of 
a dSph’s global SFR to its properties on the galactic scale. 
While we recognize the important role of molecular gas in 
star formation, our simple approach here will focus on the 
net cold gas in atomic and molecular forms, which dom¬ 
inates the total gas mass in the star-forming disk. Based 
on observations of seven nearby spiral galaxies resolved to 


750 pc, Bigiel et al. (20081 found that over the range of 
Eh ~ 1-10^ Mq pc“^ for the net surface mass density of H 
atoms and H 2 molecules in star-forming regions, Esfr can 
be described by 

1 . 85 ± 0.70 


Esfr 


Mq yr“i kpc 


= 10 


- 2 . 39 ± 0.28 


Eh 


10 Mq pC" 


.( 10 ) 


More recently, Roychowdhury et al. (20151 carried out a 


spatially-resolved study of star-forming regions dominated 
by H atoms in nearby massive spiral and faint dwarf irreg¬ 
ular galaxies. They found a much tighter relation 


Esfr 


Mq yr“i kpc 


10 ' 


10 Mq pc“ 


( 11 ) 


over the range of 'Eg ~ 0.3-30 Mq pc“ . In addition, the 
Kennicutt-Schmidt law in Eg. (J^ describes the global Esfr 
for individual galaxies with Eg ~ 10-10^ Mq pc“^. The 
above three empirical results on SFLs have overlapping 
ranges of Eg and are in broad agreement over these ranges 
within their uncertainties. Guided by these results, we as¬ 
sume that the global SFR 'i/)(t) in a dSph is related to the 
total mass Mg{t) of gas in its star-forming disk as 


, , , dM* , , , 


Mgit) 

Mq 


( 12 ) 


where A»(t) is a rate function. 

We can use the effective area Adisk(t) of the star-forming 
disk to rewrite Eg. (]12[) as 


EsFR(t) 


A*(t) 


Mq yr“i kpc 


101-®“ Mq yr- 


Adisk (t) 
10 kpc^ 


^g(t) 


10 Mq pC“ 


(13) 


For the above equation to agree with the empirical SFLs in 
Eqs. (J^, ( [To] ), and ( [11] ), we take 

0.-1 


A*(t) ~ 10' 


10 kpc 

Adisk(t) 


Mq yr" 


(14) 


and examine cases of q = 1, 1.5, and 2. We assume that 
Adisk(t) grew with time until a dSph became a satellite at 
t = tsat and then remained fixed at Adisk(tsat) during the 
subsequent SFH. As the majority of stars in the present- 
day Fornax are distributed within a region of r* ~ 2 kpc in 
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Figure 1. Data Boer et al.||2012b[ filled circles with error 
bars) on Fornax’s global star formation rate -0 as a function of 
time t since the big bang. The data are binned according to time 
(converted from stellar age) and the horizontal error bar repre¬ 
sents the size of each time bin. The filled circles give the mean 
for the average star formation rate -0 in each bin and the vertical 
error bar represents the Icr uncertainty in 0. A quadratic-spline 
fit (see Appendix 0 corresponding to the filled circles is shown 
as the solid curve. We take the centroid of the first time bin at 
t Ri 0.29 Gyr as the onset and t ^ 13.6 Gyr given by the fit as 
the end of star formation in Fornax. 


radius |de Boer et al. 2012b|, we take Adisc(isat) ~ ~ 

10 kpc^ and estimate A*(tsat) ~ ]^Q-2-8a yr“^. We will 

provide another justification for the choice of A*(faat) in 
Note that the case of q = 1 is special because the corre¬ 
sponding A* (t) is independent of time [see Eq. (141]. In this 
case, A* ~ 10“^° Mq yr”'^ is simply the fixed rate of gas 
consumption by star formation. 


3 GAS MASS, NET GAS FLOW, AND HALO 
EVOLUTION FOR FORNAX 


The data of de Boer et al. (2012b I on Fornax’s are 


binned across 13.75 Gyr and give the average SFR for 
each time bin. As discussed in Appendix we obtain a 
smooth ip{t) using a quadratic-spline fit that conserves the 
total number of stars formed in each bin and guarantees the 
continuity of ’tp{t) and dtp/dt. The fitted corresponding 
to the mean for ip in each bin is shown along with the data 
in Fig. [2 

By choosing an appropriate A,(f), we estimate Mg{t) 
from Eq. (|12|| as 


Mg{t) 

Me 


V’(^) 

A.(f) 


l/o 


(15) 


We first fix \*(t) = A*(tsat). As discussed in [2.3 compar¬ 
ison of our assumed SFL with empirical results gives an 


estimate of At(tsat) ~ 10 


’ Me yr . We then seek to 


identify fgat, the time at which Fornax became an MW satel¬ 
lite, from the evolution of the net gas flow rate 


. , dM„ ,, . 

^F{t) = ^+iP{t). 


(16) 


A net inflow onto or outflow from the star-forming disk of 



Figure 2. Total mass Mg for cold gas in Fornax’s star-forming 
disk as a function of time t derived from Eq. ID and the smooth 
0(t) shown in Fig. ^for the baseline case of a = 1.5 and A* (tsat) = 
10“^^ Mq yr“^. The top curve assumes a disk of fixed size with 
A*(t) = A* (tsat). The middle (bottom) curve is for a disk growing 
until tsat ~ 4.8 (7.8) Gyr. The middle and top (bottom) curves 
have the same Mg{t) for t > 4.8 (7.8) Gyr. Note that all curves 
reach Mg = 0 at t Rt 13.6 Gyr, the end of star formation in 
Fornax. 


Fornax corresponds to AF{t) > 0 or AF{t) < 0, respec¬ 
tively. We expect that the net gas inflow would be suppressed 
at t ~ tsat due to gas loss through ram-pressure stripping 
and tidal interaction with the MW ( Mayer et al.|[2006 see 
also [j^. 

Crucial to the identification of tsat through features of 
AF{t) are the uncertainties in ip{t). Ideally, the permitted 
smooth ip{t) should be obtained from the raw observational 
data by taking into account correlated errors in each time 
bin. Such an analysis is beyond our scope here. Upon ad¬ 
vice from T. J. L. de Boer, we give instead an approximate 
treatment as follows. We first make a Monte Carlo realiza¬ 
tion of the data by sampling ip for each time bin according to 
a Gaussian distribution with the mean and standard devia¬ 
tion specified by the corresponding mean and Ict uncertainty 
shown in Fig. The realization is accepted only if it gives 
(3.12 ± 0.05) X Me for the total mass of stars formed 
(personal communications from T. J. L. de Boer). Out of 
a total of 10^ realizations, only 3806 satisfy the constraint. 
For each accepted realization, we obtain a quadratic-spline 
fit to ip{t) and use it to calculate Mg{t) and AF{t) from 
Eqs. (151 and (161, respectively. 


3.1 Baseline case 

We first discuss our baseline case of a = 1.5 and At(taat) = 
10“^"^ Me yr~^- The Mg{t) calculated from the smooth ip{t) 
shown in Fig. assuming that A,(t) = A*(tsat) is shown as 
the top curve in Fig. The corresponding AF{t) is shown 
as the curve in Fig. |^. The shaded region in this figure 
indicates the 68% confidence interval for AF{t) estimated 
from the Monte Carlo simulations constrained by the total 
mass of stars formed. 

We take tsat to coincide with the onset of suppression 
of the net gas inflow. Based on the AF{t) shown in Fig. 
we consider two possibilities of tsat ~ 4.8 and 7.8 Gyr, re¬ 
spectively. Using the Via Lactea II cosmological simulations. 
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Figure 3. Rate AF of net gas flow to or from Fornax’s star¬ 
forming disk as a function of time t estimated from Eqs. El 
and and the data shown in Fig. for the baseline case of 
a = 1.5 and A*(isat) = 10“^^ Mq yr“ . A net inflow or outflow 
corresponds to AF > 0 and < 0, respectively. In each panel, the 
curve is calculated from the smooth ip(t) shown in Fig. and 
the shaded region indicates the 68% confidence interval. Panel 

(a) assumes a disk of fixed size with A*(t) = A*(tsat). Panels 

(b) and (c) are for a disk growing until tgat ~ 4.8 and 7.8 Gyr, 
respectively. We attribute the onset of suppression of the net gas 
inflow at t Ri 4.8 or 7.8 Gyr to effects from Fornax’s becoming an 
MW satellite. See text for details. 


Rocha et al. (20121 found that the infall times tinfaii for 


the MW satellites can be constrained by their present-day 
kinematics (radial velocities, proper motions, and Galacto- 
centric positions). They fonnd that tinfaii ~ 4.8-5.8 or 7.8- 
9.8 Gyr for Fornax, with the former range being more prob¬ 
able. It is intriguing that our estimates of fsat based on For¬ 
nax’s SFH are consistent with theirs for tinfaii based on For¬ 
nax’s present-day kinematics. 

For each choice of fsat, we assume that the Fornax halo 
stopped growing at t = fsat and that at t < tsat, its total 
mass Mh{t) followed the median halo growth history in the 
ACDM cosmology. As discussed in Appendix we use the 
model of Zhao et al. ( 2009 1 to estimate the dependence of 
the CDM density profile on the halo mass and collapse time, 
and select the appropriate Mh{tsat) to obtain a mass M{< 
r' 1 / 2 ) = 7.39lQ'3g X 10^ Mq enclosed within the half-light 


radius ri 72 = 944±53 pc as derived from observations (Wolf 


et al. 


_ ^20101. We find that Mh{tsat) ~ 1-8 x 10^ Mq and 

2.4 X 10^'Mq for tsat ~ 4.8 and 7.8 Gyr, respectively. Using 
CDM mass distributions in subhalos from high-resolution 
simulations of the Aquarius Project, Strigari et al. (2010 1 



t (Gyr) 


Figure 4. Total mass Mh of the Fornax halo as a function of time 
t before it became an MW satellite at t = teat - The solid (dashed) 
curve is for taat ~ 4.8 (7.8) Gyr. These results are estimated 
from the median halo growth history in the model of | Zhao et al.| 
( |2009[| and the present mass within the half-light radius as derived 
by [Wolf et al.| ( |2010^ from observations (see Appendix ^ . The 
dotted curve shows the characteristic halo mass Mc(f) used in 
Eq. l|18|, which gives the mean baryonic mass fraction for halos 
in a reionized universe based on simulations of lOkamoto et al.l 

< |2008l l. 


found a present-day halo mass of 7 x 10° Mq for Fornax 
based on its stellar kinematics. Considering the approximate 
nature of our approach and that the Fornax halo likely lost 
CDM at f > tsat due to tidal interaction with the MW (see 
sQ, we regard our estimates of Mh{tsat), especially the lower 
one for tsat ~ 4.8 Gyr, as consistent with this more rigorous 
result on Fornax’s present-day halo mass. 

Using the model of |Zhao et al. (20091, we show Mh{t) 
at t ± tsat for tsat ~ 4.8 and 7.8 Gyr as the solid and dashed 
curves, respectively, in Fig.|^ We now make better estimates 
of Mg{t) for t < tsat by considering a growing star-forming 
disk. In accordance with the model of |Mo, Mao fc White| 
(19981, we assume that the effective disk radius grew in pro¬ 
portion to the virial radius Tvir of the halo (see Appendix [C|) 
until t = tsat and then remained fixed during the subsequent 
SFH. Thus we obtain 

\ f A*(fsat)[rvir(Gat)/rvir(t)]^^““^^ t<tsat, 

X,{tsat), t^tsat. ^ 

For the baseline case of a = 1.5 and A*(fsat) = 

10“^"^ Mq yr“^, we show the evolution of Tvir t <C tsat 
for the two estimates of tsat in Fig. The corresponding 
Mg{t) calculated from the smooth il’{t) in Fig. [^is shown 
as the solid and dashed curves in Fig. The correspond¬ 
ing AF{t) is shown in Figs. and|^. Note that a growing 
disk does not change AF{t) at t < tsat qualitatively, and 
therefore, does not affect our above estimates of tsat- 

The choice of A* (t) can now be justified by comparing 
our results on Mg{t) and Mh{t) at t < tsat with simulations 
of gas accretion by halos after reionization of the intergalac- 
tic medium (IGM) at time treion- Okamoto, Gao & Theuns 
( |2008[ ) showed that at f freion, the baryonic mass fractions 
for halos of mass M are distributed around the mean value 


MM,t) = (/,){! ± (2^/® - l)[Ma{t)/Mf} 


- 3/2 


(18) 


where Mc(t) is a characteristic halo mass determined from 
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Figure 5. Virial temperature Tvir (solid curve) and radius rvir 
(dashed curve) for the Fornax halo as functions of time t before 
it became an MW satellite at t = fsat- Panels (a) and (b) are for 
taat ~ 4.8 and 7.8 Gyr, respectively. The horizontal part of the 
solid curve corresponds to the transition from a neutral to a fully 
ionized primordial gas at Tvij. = 1.5 X 10"* K. See Appendixfor 
details. 


their simulations and is shown as the dotted curve in Fig. 
(/i,) = r2i,/(f2cDM+flt.) ~ 0.16 is the cosmic mean mass frac¬ 
tion of baryons, and flcDM = 0.243 and fli, = 0.047 are the 
fractional contributions to the critical density of the present 
universe from CDM and baryons, respectively. The above 
equation gives fb = (fb) for M = Mc{t) and fb = (fb) 
in the limit M ^ Mc{t). Using M = MfAt) for t < 
we calculate fb{t) = fb(M,t) from Eq. (181 and show the 
result as the dotted curve in Fig.|^ This result can be com¬ 
pared with fg{t) = Mg(t) / Mh{t) for our baseline case with 
tsat ~ 4.8 (7.8) Gyr, which is shown as the solid (dashed) 
curve in Fig. &)• 

In comparing the fg{t) and fb{t) shown in Fig. we 
note that Okamoto et al.| (20081 assumed a reionization red- 
shift of « = 9 corresponding to t = 0.55 Gyr. We now esti¬ 
mate the time treion for reionization of the IGM surrounding 
Fornax. The oldest stellar age bin for the data on Fornax’s 
SFH is 13-14 Gyr. In view of the substantial uncertainties 
of ~ 2.5 Gyr for these ages ( de Boer et al.|[2012b (, we take 
the centroid of this age bin at 13.5 Gyr as a reasonable age 
estimate for Fornax’s oldest stars. With an adopted age of 
13.79 Gyr for the universe, we estimate that star formation 
started in Fornax at tsF ~ 0.29 Gyr (resdhift zsp « 14.5), 
when the halo mass reached Mhitsp) ~ 1.3 x 10® Mq or 
1.1 X 10® Mq for tsat ~ 4.8 or 7.8 Gyr, respectively (see 



t (Gyr) 


Figure 6. Fornax’s gas mass fraction fg = Mg/M^ as a function 
of time t (solid and dashed curves) before it became an MW 
satellite at t = tsat for the baseline case of a = 1.5 and A*(tsat) = 
10“^^ Mq yr“^. Panels (a) and (b) are for tsat ~ 4.8 and 7.8 Gyr, 
respectively. In each panel, the shaded region indicates the 68% 
confidence interval due to uncertainties in Mg{t) only. The dotted 
curve shows the mean baryon mass fraction fh{t) calculated from 
Eq. ( |18[ | for Fornax-like halos in a universe reionized at ^ = 9 
{t 5^ 0.55 Gyr). The horizontal dotted line indicates the cosmic 
mean mass fraction of baryons {/b}. We estimate Geion ~ 0.52 
and 0.49 Gyr corresponding to ^reion ~ 9-4 and 9.8 for panels 
(a) and (b), respectively, and set fg{t) = (fb) for t < treion- The 
agreement between fg and fb justifies the choice of A*(t). 


Fig.0. This is reasonable because the corresponding virial 
temperature T„ir{tsF) ~ 1-8 x 10® or 1.6 x 10® K (see Ap¬ 
pendix and Fig. could readily enable cooling by H 2 
molecules as required for formation of the first stars (see e.g., 
Bromm fc Larson|2004| for a review). With tsF ~ 0.29 Gyr, 
star formation most likely started in Fornax at t < treion. For 
tsF < t < treion, accretion of gas by the Fornax halo should 
be proportional to that of CDM and fg{t) = Mg{t)/Mi,{t) 
should be close to the cosmic mean mass fraction of baryons 
(fb) « 0.16. The fg{t) shown in Fig. ([^) reaches (fb) at 
t ~ 0.52 (0.49) Gyr and exceeds this value at earlier times 
(not shown). Our assumed SFL most likely does not apply to 
the early epoch near the onset of star formation. So we could 
take Ueion ~ 0.52 (0.49) Gyr corresponding to Zreion ~ 9.4 
(9.8) and make a revised estimate of Mg{t) ~ {fb}Mh{t) for 
tsF < t < treion- The above estimates of Zreion are in good 
agreement with « = 8)(® found for reionization of the IGM 
surrounding MW dwarf galaxies in simulations by|Busha et| 

poTo) ). 

With Zreion ~ 9.4 (9.8) being close to the value as¬ 
sumed by Okamoto et al. (20081, our estimated fg can be 
directly compared with fb and the agreement is rather good 


al. 
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Figure 7. Same as Fig. Hji but for the case of a = 1 and 
A*(£sat) = 7 X 10“^^ M0 yr“^. Note that A*(t) = A*(tsat) is 
independent of time for this case. 



t (Gyr) 


Figure 8. Same as Fig. but for the case of a = 1 and 
A*(tsat) = 7 X 10“^^ Mq yr“i. Note that A*(f) = A*(tsat) is 
independent of time for this case. 


as shown in Fig. We note that this agreement is not af¬ 
fected by adding the contribution of stars to fg as the mass 
of stars is always significantly smaller than that of gas for 
t < tsat- Therefore, we consider that our choice of A*(t) with 
A*(tsat) = 10“^“^ Mq yr“^ for the baseline case is reasonable. 


3.2 Other cases 


The case of a = 1 is special because A*(t) = A*(tsat) is in¬ 
dependent of time (see (2.31. We choose A*(tsat) for this 
case to obtain the same Mg at t ~ 4.8 Gyr as the top 
curve in Fig. for the baseline case. This gives A*(tsat) = 
7 X Mq yr”'^, which is consistent with the estimate 

based on empirical SFLs as discussed in §2.3| The corre¬ 
sponding AF{t) is shown in Fig. where the same on¬ 
set of suppression of the net gas inflow at faat ~ 4.8 or 
7.8 Gyr as for the baseline case can be seen. In Fig. we 
take tsat ~ 4.8 Gyr and show the corresponding fg{t), which 
gives treion ~ 0.88 Gyr (^reion ~ 6.3). As fb assumes a = 9 for 
reionization, straightforward comparison between fg and ft 
can be made only for t > 0.88 Gyr. The agreement is rather 
good for t > 1.2 Gyr. 

For the case of q = 2 , we first consider A*(t) = A*(taat) 
and choose A*(taat) to obtain the same Mg a.t t ~ 4.8 Gyr 
as the top curve in Fig. for the baseline case. This gives 
A*(tsat) = 1.4 X 10~^® Mq yr“^ in accord with the esti¬ 


mate based on empirical SFLs (see (2.31. The corresponding 
AF{t) shown in Fig.[^ again suggests faat ~ 4.8 or 7.8 Gyr. 
The AF{t) for a disk growing until t « 4.8 Gyr is shown in 


Fig.j^ and the corresponding fg{t) in Fig. 10 An estimate 
of treion ~ 0.36 Gyr (zreion ~ 12.3) Can be made and there 
is good agreement between fg and fb at t > 0.55 Gyr, i.e., 
after the time of reionization assumed for the latter. 

Based on the above results, we consider that our results 
on faat, aud hence Mb.{tsa,t), are insensitive to the exact form 
of our assumed SFL so long as it is consistent with the em¬ 
pirical SFLs discussed in §2.3| 



Figure 9. Same as Figs.|^ and[^ but for the case of o = 2 and 
A*(tsat) = 1.4 X 10“^® Mq yr“i. 


4 GLOBAL GAS DYNAMICS AND 

EVOLUTION OF FORNAX AS AN MW 
SATELLITE 

Based on the discussions in ii the baseline case appears to 
give representative results on the general evolution of For¬ 
nax. Below we focus on this case with fsat ~ 4.8 Gyr to 
discuss the global gas dynamics and evolution of Fornax as 
an MW Satellite. Considering that our discussion is mostly 
qualitative, we use the solid curve for AF{f) at t > 4.8 Gyr 
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t (Gyr) 


Figure 10. Same as Fig. but for the case of o = 2 and 
A«(tsat) = 1-4 X 10“^® Mq yr“i. 

in Fig. Ep as a guide. This is reproduced in Figure pT] :, which 
shows that AF dropped to zero at t ~ 6 Gyr and the first 
net outflow {AF < 0) occurred for t ~ 6-7 Gyr. This was 
followed by the final round of net inflow {AF > 0) for t « 7- 
8.6 Gyr, and then by several episodes of net outflow until 
star formation ended in Fornax. As we discuss below, the 
episodic gas flows at t > 4.8 Gyr shown In Fig. [nj; were 
driven by Fornax’s orbital motion and its tidal interaction 
with the MW. 


4.1 Orbital motion and tidal interaction 


Observations show that Fornax currently moves in an ellip¬ 
tical orbit with a period of ~ 3.2 Gyr. Its distance from 


Ra ~ 152 kpc at the apogalacticon i Piatek et al.|2007 l. Us¬ 
ing the median halo growth history and the corresponding 


halo structure given by the model of Zhao et al. (20091, we 


find that for a present-day MW halo mass in the range of 
(1-2) X10^^ Mq, the MW mass enclosed by Fornax’s present- 
day orbit is ~ 1.5 times the MW halo mass at taut ~ 4.8 Gyr 
when Fornax became a satellite. In view of this moderate 
change, we assume that Fornax’s orbit for t > 4.8 Gyr can be 
approximately described by the present-day parameters. For 
estimates, we take the MW mass enclosed within Fornax’s 
orbit to be ~ 7.5 x 10^^ Mq on average. This gives an 

average orbital velocity of upor ~ {GM^^/a)^^^ ~ 150km/s 
for Fornax, where a = {Rp + Ra)/2 ~ 135 kpc. For compari¬ 
son, the circular velocity characterizing matter motion inside 
the Fornax halo at fsat ~ 4.8 Gyr is Udrc = (GM^/rvir)^^^ ~ 
21 km/s with Mh ~ 1.8 x 10® Mq and Pvir ~ 17 kpc (see 
Fig.§. 

The effects of tidal interaction between the MW and 
Fornax can be estimated from the radii rns.p and rns.a of the 
instantaneous Roche spheres | |King|[T962[ ) at the perigalac- 
ticon and apogalacticon, respectively. Ignoring the small ec¬ 
centricity of ~ 0.13 for Fornax’s orbit ( Piatek et al.|[2007 l, 
we have 

)"RS,P ^ r-HS,a ^ ( Mpor \ ^ 


Rv 


Ra 




(19) 


where Mpor is the total mass of Fornax and generally differs 
from Mh at tsat ~ 4.8 Gyr due to tidal interaction. The CDM 



Figure 11. Evolution of the Fornax halo as an MW satellite 
at i > 4.8 Gyr. Its total mass Mpor (a) and Roche-sphere radius 
rjTs (b) are estimated at the apogalacticon (indicated by “a”) and 
perigalacticon (indicated by “p”), and connected by line segments 
to guide the eye. For convenience of discussing the effects of tidal 
interaction, the relevant AF{t) for the baseline case (c) is also 
shown. See text for details. 


and gas outside the Roche sphere would be lost from Fornax. 
However, because the Roche sphere would expand as Fornax 
moved from the perigalacticon to the apogalacticon, some of 
the GDM and gas lost earlier could be reaccreted. Guided by 
the simulations of Han et al. ( 2012 I and our derived AF, we 
assume that reaccretion occurred only once. Subsequently, 
Mpor settled down to a value set by the mass enclosed within 
fRS.p- For estimates, we use the density profile given by the 
model of Zhao et al. (20091 for a halo with Mh ~ 1.8 x 
10® Mg at taat ~ 4.8 Gyr. Assuming that Fornax was at the 
apogalacticon at taat ~ 4.8 Gyr with Mpor ~ 1.8 x 10® Mq, 
we estimate that its total mass was quickly reduced to ~ 
1.6 X 10® Mq , which corresponds to the mass enclosed within 
the Roche sphere at this time. The mass of the Fornax halo 
was further reduced to ~ 1.3 x 10® Mq at t ~ 6.4 Gyr when 
it was at the perigalacticon with a smaller Roche sphere. It 
then reaccreted some of the CDM (and gas) lost earlier to 
obtain Mpor ~ 1.5 x 10® Mq at t ~ 8 Gyr when it was at 
the apogalacticon again. After that Mpor settled down to ~ 
1.2 X 10® Mq, which is enclosed within trs.p estimated from 
Eq. (191 for this value of Mpor- During the above evolution 
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of the Fornax halo, rns,a started at ~ 14 kpc and settled 
down to ~ 12 kpc while rns.p decreased from ~ 11 kpc to 
~ 10 kpc. The evolution of Mpor, rRs.a, and rpis,p is sketched 
in Fig. 

Based on the above discussion, we estimate that Fornax 
lost ~ 1/3 of its CDM through tidal interaction. However, 
this loss occurred in the ontermost region of the halo and 
had little effect on the strnctnre of its interior. For example, 
because ri /2 ~ 944 pc was well within the Roche sphere at 
all times, the mass M(< ri/ 2 ) enclosed within ri/ 2 , which is 
used to determine Mh at tsat ~ 4.8 Gyr (see Appendix m, 
would not have been affected by tidal interaction. Similarly, 
with a radius of ~ 2 kpc, the star-forming disk would have 
remained largely intact as Fornax orbited the MW so long 
as there was a sufficient supply of cold gas. 


4.2 Gas inflow and outflow 


Over Fornax’s SFH, the star-forming disk was surrounded 
by hot gas either accreted from the IGM or expelled from the 
disk due to feedback from star formation, which includes SN 
explosions and heating by stellar radiation. Cooling of this 
hot gas gave rise to an inflow onto the disk while stellar feed¬ 
back produced an outflow from it. The net gas flow rate AF 
represents the difference between the two. At t < 4.8 Gyr, 
accretion of the IGM by the growing Fornax halo dominated 
its global gas dynamics and supplied hot gas that cooled to 
provide a net inflow. In general, reducing the amount of 
hot gas would suppress the inflow, thereby causing AF to 
decrease. This accounts for the decline of AF for t « 4.8- 


6 Gyr (see Fig. Ill as accretion of the IGM was disrupted 


and the initial tidal interaction between Fornax and the MW 
resulted in loss of hot gas along with CDM. This gas loss 
was also aided by ram-pressure stripping as discussed in 
§4.3| The corresponding suppression of the inflow allowed 
the brief hrst occurrence of a net outflow for t ~ 6-7 Gyr 
as hot gas escaped from the star-forming disk. This outflow 
carried only ~ 2 x 10® Mq and was soon quenched when 
~ 10^ Mq of the hot gas lost earlier was reaccreted and 
cooled to produce an inflow for t « 7-8.6 Gyr (see Fig. 111. 


As mentioned in §4.1[ some of the gas and CDM lost 
earlier was reaccreted because the Roche sphere expanded 


at the apogalacticon (see Fig. 111. This reaccretion was dis¬ 
cussed for gas by Nichols, Lin &: Bland-Hawthorn|( 20121 and 


seen for CDM in simulations by Han et al. ( 2012[). However, 


our derived AF indicates that there was only one episode 
of signihcant reaccretion for Fornax, which might represent 
a variation from the scenario of episodic starbursts with re¬ 
peated gas expulsion and reaccretion as discussed by |NichoIs| 
et al. (20121. In the absence of further reaccretion, several 
episodes of net outflow occurred in Fornax for t « 8.6- 

13.6 Gyr (see Fig. |ll[ |. We consider that the episodic nature 
of these outflows is mainly due to Fornax’s orbital motion 
and is very different from the episodic behavior purely driven 
by stellar feedback as discussed for example, by [Stinson et| 
al. (20071. Specifically, we argue that although stellar feed¬ 


back initiated the escape of hot gas from the star-forming 
disk, the loss of such gas from the Fornax halo at t « 8.6- 

13.6 Gyr was driven by tidal interaction and ram-pressure 
stripping. 


4.3 Ram-pressure stripping 

We consider that the gas subject to ram-pressure strip¬ 
ping was dispersed outside Fornax’s star-forming disk but 
within its Roche sphe re of ~ 10 kpc in radius. From the 
we estimate that ~ 10^ Mq of gas 


11 


AF shown in Fig 
were stripped for each episode. So this gas had a density of 
n ~ 10“^ cm“® while inside Fornax. Its internal pressure 
was P ~ nkT ~ rupriv^i-^, with T being its temperature, 
Uth ~ {kTjmpY^^ ~ 29(r/10® km/s its thermal ve¬ 

locity, and k the Boltzmann constant. The ram pressure can 
be estimated as 


P 2 

4 ram ^^p^med^For? 


( 20 ) 


where rimed is the density of the gaseous medium through 
which Fornax moved. Ram-pressure stripping occurred when 
Dram exceeded P, i.e.. 


2 2 
^med^For ^ ^^th* 


( 21 ) 


For the hot gas accreted from the reionized IGM right be¬ 
fore Fornax became an MW satellite, we estimate Uth ~ 
Wcirc ~ 21 km/s. This gas would be stripped for rimed ^ 
2 X 10"® cm"®. 

The majority of the gas stripped from Fornax first es¬ 
caped from its star-forming disk due to stellar feedback. 
For a crude estimate, we assume that the gas in this disk 


consisted of three components (McKee & Ostriker 19771: 
cold gas with a density of ricoid ~ 40 cm"'” and a tem¬ 
perature of Tcoid ~ 80 K, warm partially-ionized gas with 


0.3 cm ® and Tw 


8000 K, and hot ionized 


gas with rihot ~ 4 x 10"® cm"® and Thot ~ 5 x 10® K. 
The total gas mass was dominated by the cold compo¬ 
nent that fed star formation. Cold gas was also converted 
by stellar feedback to maintain the hot component as hot 
gas escaped from the disk in an outflow. Assuming adia¬ 
batic cooling during its escape, we estimate an upper limit 
of T ~ 2hot(n/nhot)®^® ~ 4 x lO’^ K corresponding to 
Uth ~ 18 km/s for the hot gas when it was dispersed inside 
the Roche sphere. Its stripping requires rimed ^ 10"® cm"®. 

Based on the above discussion, with rimed ^ 2 x 


bit at ~ 118-152 kpc from the MW center, ram-pressure 
stripping would have occurred to both accreted hot gas for 
t « 4.8-6 Gyr and gas heated by stellar feedback for t ~ 8.6- 

13.6 Gyr in Fornax. The gas density in the present MW halo 
can only be estimated by indirect means and its time evo¬ 
lution over the MW’s history can only be estimated by sim¬ 
ulations of galaxy formation. Simulations of jKaufmann et| 

at ~ 100-200 kpc 


al. 


(20091 give rimed ~ 10 


from the galactic center with no significant changes over 
a period of 10 Gyr. These results are consistent with esti¬ 


mates based on observations of MW halo clouds (e.g., Hsu 
jet al. ]20TT| . However, [Muralij ( |2000| derived a present limit 
of rimed < 10"® cm"® at 50 kpc from the MW center by 
considering effects on the Magellanic Stream. For compar¬ 
ison, the cosmic mean baryon density at tsa,t ~ 4.8 Gyr 
(2sat ~ 1.33) was {rib) « 3 x 10"® cm"®. We consider 
rimed ~ 3 X 10"®-10"® cm"® as reasonable around Fornax’s 
orbit for t > 4.8 Gyr. Therefore, ram-pressure stripping 
would have been effective in removing gas from Fornax. 

For the three gas components in Fornax’s star-forming 
disk with densities and temperatures assumed above, they 
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were in pressure equilibrium and their internal pressure was 
> 10^ times the ram pressure. Consequently, for those con¬ 
ditions the ram pressure would have had little impact on the 
disk. However, when a sufficient amount of gas had been re¬ 
moved, the three-component structure would be disrupted 
and gas density in the disk would drop, eventually allowing 
gas to be ram-pressure stripped directly from the disk. This 
might account for the sharp increase in the net outflow rate 
at the end {t ~ 12.8-13.6 Gyr, see Fig |ll| :). 


4.4 Gas loss from Fornax 


In summary, we propose the following outline for gas dy¬ 
namics in Fornax as it orbited the MW. Tidal interaction 
was responsible for removing ~ 1/3 of its original CDM. 
Because collisionless CDM does not suffer ram pressure, the 
CDM lost during Fornax’s first orbital period appeared to 
have stayed close to the orbit and kept the simultaneously- 
lost gas from dispersing. This was crucial to the reaccretion 
of some of the lost CDM and gas as the Roche sphere ex¬ 
panded near the end of this period. Subsequently, there was 
no significant clustering of CDM near the orbit and once 
ram-pressure-stripped gas moved outside the largest Roche 
sphere of ~ 12 kpc in radius, it was lost from Fornax. We 
note that although the change in the Roche-sphere radius 
from ~ 10 kpc at the perigalacticon to ~ 12 kpc at the 
apogalacticon appeared small, the corresponding change in 
the Roche-sphere volume was by a significant factor of ~ 1.7. 
The non-monotonic evolution of net gas outflow for t « 8.6- 
13.6 Gyr might have resulted because an expanding Roche 
sphere would oppose while a shrinking one would enhance 
the effect of ram pressure. This would explain the coinci¬ 
dence of AF ~0att~11.2 Gyr with the transition from 
an expanding Roche sphere to a shrinking one (Fig. |ll[ ). We 
also note that subsequently the net outflow rate lAFI in¬ 
creased as Fornax approached its perigalacticon. However, 
the net outflow rate started to decrease at t ~ 12.3 Gyr 
and dropped to ~ 0 when Fornax reached its perigalacticon 
at t ~ 12.8 Gyr. This is in contradiction to the expected 
effect of the smallest Roche sphere at the perigalacticon. In¬ 
terestingly, several observations (e.g., Goleman et al.|[2004 


2005 de Boer et al.|2013 1 found evidence that Fornax might 
have accreted a smaller satellite or reaccreted some of the 
lost gas ~ 1.5-2 Gyr ago {t ~ 11.8-12.3 Gyr). This would 
have explained the puzzling drop of the net outflow rate de¬ 
scribed above. In any case, it appears that Fornax lost all of 
its gas at t « 13.6 Gyr before returning to its apogalacticon 
(estimated to occur at t ~ 14.4 Gyr or ~ 0.6 Gyr from now). 

Gas loss from satellites through ram-pressure stripping 
and tidal interaction was studied in detail by simulations of 


Mayer et al. (20061, who emphasized the importance of the 


UV background in keeping the gas widely distributed inside 
the satellites. Gas depletion in Local Group dwarfs was stud¬ 


ied analytically by Nichols & Bland-Hawthorn (20111, who 


emphasized heating of gas by stellar feedback. In general ac¬ 
cord with these detailed studies, we propose that gas in the 
star-forming disk was first heated by SNe to Thot ~ 5x 10® K, 
and then escaped to be dispersed outside the disk before get¬ 
ting removed from Fornax. It is also possible that heating 
by both SNe and the UV background was needed, but we 
focus on SN heating below. 

Using the same initial mass function as adopted by 


de Boer et al. 1 2012b| in deriving the SFR and assuming 


that stars of 8-120 Mq result in SNe, we obtain an SN rate 
of Rs^it) ~ 10~'^Ip{t)/MQ. For SNe to heat a total mass 
AMhot of gas that escaped during an episode of outflow, the 
required efficiency per SN is 


(AMhot/mp)fcrhot 

AVsN^^expl 



where AAsn is the total number of SNe occurring in this 
episode and Fexpi is the average explosion energy per SN. 
From the and AF shown in Figs, [^and EH respectively, 
we estimate that AMhot ~ 1.2 x 10^ Mq (2.2 x 10^ Mq) and 
ANsn « 6 X lO’^ (1.4 X 10"*) for the episode of t « 8.6- 
11.2 Gyr (11.2-12.8 Gyr), which correspond to esN ~ 1.7% 
(13%) for Thot ~ 5x 10® K and Fexpi ~ 10®'^ ergs. For the last 
episode of t ~ 12.8-13.6 Gyr, there was a sharp increase in 
the net outflow rate and AMhot ~ 8 x 10® Mq, which would 
have required esN ~ 51% with ANsn ~ 1-3 x 10®. As noted 
in §4.3| for this terminal episode, gas could have been ram- 
pressure stripped directly from the star-forming disk due 
to the decrease in density. So a lower SN-heating efficiency 
could have been sufficient. Heating by Type la SNe, which 
is ignored in the above discussion, could also reduce the re¬ 
quired esN- Therefore, it appears that SN-heating coupled 
with ram-pressure stripping and tidal interaction for a few 
orbital periods was responsible for removing gas and termi¬ 
nating star formation in Fornax. 


5 DISCUSSION AND CONCLUSIONS 


We have presented an empirical model for the halo evolu¬ 
tion and global gas dynamics of Fornax. Guided by data 
on star formation, especially those reported by [Bigiel et al.| 
(20081 and Roychowdhury et al. ( 2015| ) for nearby galaxies, 
we have assumed a global SFR ipit) = A*(t)[Mg(t)/MQ]“ 
that is dependent on the total mass Mg{t) for cold gas in the 
star-forming disk. We have examined three different cases 
of q: = 1, 1.5, and 2, and chosen the corresponding A*(t) 
in agreement with the empirical SFLs. For each case, we 


have used the data on Fornax’s i/>(t) provided by de Boer 


et al. (2012bI to derive Mg{t) and estimate the rate of net 


gas flow AF{t) to or from Fornax’s star-forming disk. We 
have identified the onset of the transition in AF{t) from a 
net inflow to a net outflow as the time tsat at which Fornax 
ceased evolving independently and became an MW satellite. 
We have determined the mass Mh{tsa,t) of the Fornax halo 
at this time by using the median halo growth history in the 
model of |Zhao et al.| ( |2009[ ) and requiring that the corre¬ 
sponding density profile give the mass enclosed within the 
half-light radius as derived from observations. We have fur¬ 
ther justified our chosen A*(t) by comparing the gas mass 
fraction fg{t) = Mg{t)/Mh{t) at t < fsat with the results 


from cosmological simulations of Okamoto et al. (2008 \ on 


gas accretion by halos in a reionized universe. 

Our main results for Fornax are the evolution of its total 
gas mass Mg{t) and net gas flow rate AF{t), the time tsat 
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when it became an MW satellite, and its halo mass M/i(tsat) 
at this time. For our baseline case oi a = 1.5 and A«(tsat) = 
10“^^ Mq yr“^, we have obtained tsat ~ 4.8 (7.8) Gyr and 
Mhitsat) ~ 1.8 X 10® Mq (2.4 X 10® M 0 ). These results are 
in broad agreement with previous studies based on Fornax’s 
orbital motion ([Rocha et al.|[2012p and stellar kinematics 
(Strigari et al. 20101. We also found that the evolution of 
the gas mass fraction fg{t) = Mg{t)/Mh{t) for t < taat is 
consistent with cosmological simulations of jOkamoto et aT] 
(20081. We have checked that these results, especially taat 
and Mfi(taat), are not sensitive to different choices of a = 1 
and 2 . 

Using the results for the baseline case with taat ~ 
4.8 Gyr and the present orbital parameters of Fornax, we 
have related AF{t) at t > 4.8 Gyr to its orbital motion as 
an MW satellite and estimated the effects of ram pressure 
and tidal interaction with the MW. The Fornax halo lost 
~ 1/3 of its GDM through tidal interaction but reaccreted 
some of the lost GDM near the end of its first orbital period. 
This reaccretion was responsible for the last episode of sig- 
nihcant net gas inflow. Otherwise, gas was removed from the 
star-forming disk as hot gas created by SN heating escaped 
and was then lost from Fornax through ram-pressure strip¬ 
ping and tidal interaction. This lasted a few orbital periods 
and eventually terminated star formation in Fornax. 

Our assumed global SFL most likely does not apply 
to the epoch near the onset of star formation and prior 
to reionization. By using the median halo growth history 
Mh(t) prior to the time treion for reionization and assuming 
the corresponding gas mass fraction fg{t) = Mg{t)/Mhit) 
to be the cosmic mean mass fraction of baryons (/b), we 
have estimated treion as the time at which the fg(t) de¬ 
rived from the assumed SFL becomes equal to (/b). We 
have obtained Leion ~ 0.52 (0.49) Gyr corresponding to 
2 reion ~ 9.4 (9.8) for the baseline case with fsat ~ 4.8 
(7.8) Gyr, Leion ~ 0.88 Gyr (zreion ~ 6.3) for the case of 
a = 1 and Geion ~ 0.36 Gyr (^reion ~ 12.3) for the case of 
a = 2. These estimates of Zreion are close to z = 8 I 2 found 
for reionization of the IGM surrounding MW dwarf galaxies 
in simulations by Busha et al. (20101. Our estimate of the 
reionization time can be improved significantly by using de¬ 
tailed CDM simulations to estimate halo evolution instead 
of the model of [Zhao et al.| ( [^009| ), which only gives the me¬ 
dian halo growth history. We note that data on SFH are 
available for other dwarf galaxies of the Local Group (e.g., 
jde Boer et al.|[2012a| |2014| jWeisz et al.||201^ . Gombining 
CDM simulations with our approach to study these systems 
may provide an interesting probe of reionization. 

In conclusion, we have provided a reasonable picture 
for the halo evolution and global gas dynamics of Fornax 
based on empirical data. The input to our model includes 
the SFH, the mass enclosed within the half-light radius, and 
the present-day orbit of Fornax. Essential to our approach is 
the SFH and we have given only an approximate treatment 
of the correlated errors in each time bin to obtain a smooth 
SFH. We emphasize that a rigorous treatment should be 
carried out to make our approach more quantitative. Our 
approach can be extended to other dSphs for which the data 
listed above are available, e.g.. Sculptor and Carina. Our re¬ 
sults on taat and Mh(taat) niay help identify the halos that 
are most likely to host the observed dSphs in CDM simula¬ 
tions of the formation of an MW-like galaxy. Our results on 


Mg{t) and AF{t) may be used along with the data on SFHs 


to develop models of chemical evolution for dSphs (e.g., Qian 
jfc Wasserburg|[2012| and references therein). Our empirical 
approach also serves to illustrate how different sets of obser¬ 
vational data may be integrated to provide a coherent de¬ 
scription of the evolution of individual dSphs. We hope that 
this would help motivate dedicated cosmological simulations 
with gas physics to reproduce the observed dSphs with de¬ 
tailed star formation and chemical evolution histories, as 
well as present-day mass distributions and orbits. Such sim¬ 
ulations are the best way to resolve the TBTF problem. 
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APPENDIX A: FITTING FORNAX’S SFR TO 
DATA 


The data of de Boer et al. (2012b I on Fornax’s SFH were 
binned according to the age t* of stars. We label the bin 
boundaries as to < H < • • • < tn in terms of the time of star 
formation t = — t*, where is the age of the universe and 

t = 0 corresponds to the big bang. We obtain a smooth SFR 
from a quadratic-spline fit to the data that conserves 
the total number of stars formed in each bin and guarantees 
the continuity of ip{t) and ijj'it) = dip/dt (see e.g., G. Ghan 
et al., in preparation). Specifically, for the ith bin 
with an average SFR ipi given by the data, we seek 


+pi{t- ti-i) + qi{t- ti-if 


(Al) 


where pi and qi are coefficients to be determined. The con¬ 
servation of the total number of stars formed in this bin 
requires 


- 1 1 2 
=i>i- -piAi - -qiAi 


where Ai = ti — ti-\. In addition, Eq. (Al I gives 

= Pi, 

'tp'di) = pi-\-2qiAi. 


(A2) 


(A3) 

(A4) 


Inspection of Eqs. (A1|-(A4| shows that 'i/)(t) for all t in 
the ith bin can be obtained once and are 

determined. 


Equation (All gives 


ipiu) = Ipiti-i) + piAi + qiAi 


(A5) 


Substituting ipiti-i) from Eq. (A21 into the above equation, 
we obtain 

ip{ti) = ipi + ^PiAi + (A6) 

All of the equations for the ith bin can be generalized to 


other bins. For example, Eq. (A2 I can be rewritten as 


- 1 1 2 
tpiU) = V’i+i - ^Pi+i^i+i - 


(A7) 


Combining Eqs. (A 61 and (A71 gives 


+ Ai J- ipi+iA,+i J- igi+i Ai+i = V’i+i - V'i-(A8) 


rPiA, 


Using Eqs. ( |A3[ ) and ( |A4[ ) to eliminate pi, qi, Pi+i, and Qi+i 
from the above equation, we obtain 


Ai ,/. . Ai 

—pj (ii_i)H- 


+ A, 


i+i, 


V’'(L)H — ’—ip'{ti+i) = ipi+i-tpi, (A9) 


which can be rewritten as a matrix equation for {'ip'{ti),0 sC 
i Si n}. The unknown {ip'{ti),0 < i < n} can be solved from 
this equation by specifying tp'ito) and xj)'{t„). 

We take xp'{to) = = 0 to obtain the smooth xp(t) 

that is shown in Fig. along with the data. We use tu = 
13.79 Gyr for the adopted cosmology to convert the stellar 
age tt into the time t of star formation. The youngest stars 
covered by the data have t* = 0.25-0.5 Gyr corresponding 
to t = 13.29-13.54 Gyr. Considering that there is no current 
star formation in Fornax, we have included 13.54 ^ t ^ 
13.79 Gyr in the fit as the nth bin with xpri = 0. To ensure 
that no net stars were formed in this bin, the corresponding 
'!/j(t) must drop from positive to negative values. We take the 
part of the fit with •ip{t) ^ 0 and ignore the unphysical part 
with xp(t) < 0. This gives an estimate oi t ~ 13.6 Gyr for 
the end of star formation in Fornax. The oldest stars covered 
by the data have = 13-14 Gyr. So formally our first bin 
corresponds to —0.21 ^ t ^ 0.79 Gyr. However, the stellar 
ages in this bin have uncertainties of ~ 2.5 Gyr (de Boer 
et al.|[^12b l. Therefore, we take the centroid of this bin at 


t = 0.29 Gyr as approximately the onset of star formation 
in Fornax and ignore the part of the fit for t < 0.29 Gyr. 
The uncertainties associated with the first and last bins for 
our fit represent end effects that are mostly restricted to the 
corresponding periods. Our main results are based on the 
fitted outside these periods. 
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APPENDIX B: DETERMINATION OF 
FORNAX’S HALO MASS 

In discussing halo evolution, it is convenient to use redshift 
2 . For a flat ACDM cosmology, we have 


dz .1 , ^ 

where 


SirGpcriz) 


Pcr{z) = Pcr(O) + z)^ + Ha] 


(Bl) 


(B2) 


is the critical density at redshift z, Pcr(O) = 2>HQ/{di'KG) 
is the critical density at the present time, Hq = 
lOO/i km s”'^ Mpc“^ is the Hubble parameter, and Q,m = 
HcDM + Hf) and Ha are the fractional contributions to Pcr(O) 
from non-relativistic matter (CDM plus baryons) and the 
cosmological constant, respectively. Throughout this paper, 
we adopt h = 0.69, Hm = 0.29 (Hcdm = 0.243, Ht. = 0.047), 
and Ha = 0.71. In calculating the growth of a halo, we use a 
primordial power spectrum with a spectral index Us = 0.96 


and the transfer function of Eisenstein & Hu (19981 with 


a present temperature of 2.726 K for the cosmic microwave 
background to obtain the linear power spectrum, the ampli¬ 
tude of which is fixed by erg = 0.82. The above cosmologi¬ 
cal parameters are consistent with the final analysis of the 
WMAP experiment ( Hinshaw et al.|2013 |. 

In addition to the total mass of a halo, another impor¬ 
tant parameter is its virial radius. For a halo of mass Mh 
collapsing at redshift z, its virial radius rvir (Bryan & Nor¬ 
man | |1998[ ) is defined through 


47r Q 

Mh = -^rvirpcr(-2)Ac, 

where 

Ac = -f 82(Hc - 1) - 39(H, - 1)^], 


and 


H, = 


Hm (1 + z)^ 


(B3) 


(B4) 


(B5) 


Hm(l + z)^ -l- Ha 

The mass distribution of a halo can be approximated by the 
Navarro-Frenk-White (NFW) density profile 

Po 


p{r) = 


(B 6 ) 


{r/rs)[l + {r/rs)f' 

where po and Ts are two parameters characteristic of the 
halo. Using the above profile, we obtain 


M. = 


AtvE p{r)dr = 47 rpor^ 


ln(l -I- c) -I- 


1 -I- c 


(B7) 


where c = Vvii/Vs is referred to as the concentration pa¬ 
rameter. For a halo of mass Mh collapsing at redshift z, the 
concentration parameter c{Mh, z) can be estimated from the 
model of Zhao et al. (20091. 


We assume that the density profile of the Fornax halo 
was approximately fixed at fgat ~ 4.8 Gyr (zaat ~ 1.33) when 
the halo mass reached Mh{zaa,t)- As discussed in Q tidal 
interaction with the MW would not have affected the density 
profile within the half-light radius ri /2 for the present-day 
Fornax. Using the NFW profile and Eq. (B71, we obtain the 
mass M{< ri/ 2 ) enclosed within ri /2 as 

M{< ri/ 2 ) = , (B 8 ) 


where 

f{x) = ln(l -I- x) -I- 


1 -(- x’ 


(B9) 


and both c and Xvir are evaluated for Mh(zsat)- Estimating 
c{Mh, Zaax) with the model of [Zhao et al. ( |2009 1, we solve 
Eq. (B 8 I by iteration to obtain Mh{zsa.t) ~ 1.8 x 10® Mq for 
Mi< ri/ 2 ) = 7.39tn:gfi X 1 O’’M 0 and ri /2 = 944 ± 53 pc 
( Wolf et al.[Mot . 


APPENDIX C: CHARACTERISTIC 
PARAMETERS OF A HALO 


We define some useful quantities (e.g., Barkana fc Loeb 
|200ip for discussing the evolution of a halo. For a halo of 
total mass Mh collapsing at redshift z, its virial radius r. 


is defined by Eq. (B31 and can be evaluated as 
1/3 


7.85 f Mh 


1 + z \^ 1 O 8 M 0 
Its circular velocity is 


IStt" 


(nil) '“’'I'”’ 


'^circ — 


GMh 


= 20.7 


Mh 


Xvir \ 10® Mq 

and its virial temperature is 


1/2 


^ kms“\(C2) 

V Xvir / 


Tvir = 


prUpV^ 


2k 

= 3.32 X 10^(1 -ba)p 

2 \ 1/3 


Mh 

108 M 0 


2/3 


(illi) 


1/3 


X 




K, 


(C3) 


where p is the mean molecular weight for those electrons, 
nuclei, and atoms that contribute to the gas pressure, and 
k is the Boltzmann constant. We take the primordial mass 
fractions of protons and '^He nuclei to be 0.75 and 0.25, 
respectively, and use p — 1.23 and 0.59 for neutral and fully 


ionized gas, respectively. Following Busha et al. (20101, we 
take p = 1.23 for Tvir < 1.5 x 10 K and p = 0.59 for 


Mh{t) for the Fornax halo shown in Fig. we show the 
corresponding rvir(t) and rvir(t) in Fig. 


© 2015 RAS, MNRAS 000,[iHT3l 















































